BRX TH-570 



Multipole Radiation in Lorentz Gauge 

Howard J. Schnitzer 1 
Martin Fisher School of Physics 
Brandeis University 
Waltham, MA 02454 



Abstract: The multipole expansion for electromagnetic radiation, valid for all wave-lengths and all distances 
from bounded sources, is presented in Lorentz gauge, rather than the usual Coulomb gauge. This gauge 
is likely to be preferred in applications where one wishes to maintain manifest Lorentz invariance. The 
presentation also serves as a useful exercise in the use of vector spherical harmonics. 

The multipole expansion is a standard issue for the description of electromagnetic radiation. 
In many applications one may employ suitable simplifications such as the long- wave limit, or re- 
striction to just the asymptotic behavior of the fields. It is of interest in other contexts to describe 
the multipole expansion, valid for all wave-lengths and all distances from bounded sources. Typ- 
ically this expansion is presented in Coulomb gauge. However for applications where one wishes 
to maintain manifest Lorentz invariance, the Lorentz gauge is preferable. Of course, the resulting 
electric and magnetic fields do not depend on gauge choice. But in some applications, the multipole 
expansion in terms of Lorentz covariant potentials is particularly useful. 

It is the purpose of these notes to describe the multipole expansion of electromagnetic ra- 
diation in Lorentz gauge valid for all wave-lengths and all distances from bounded sources, as we 
are not aware that this is available elsewhere. Other useful discussions, with other aspects of the 
subject, are to be found in refs. [1—5]. Various useful properties of spherical Bessel functions are 
available in standard references on quantum mechanics, for example. 

The electric and magnetic fields may be described in terms of vector and scalar potentials A 
and (p, 

H(r, t) = V X A(r, t) (1) 



and 

at 

In Lorentz gauge, one has the gauge condition 



E(r, t) = -V0(r, A(f, t) . (2) 



V-l+if = 0, (3) 
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which implies the potentials satisfy the wave-equations 



- V 2 A(r, t) + £ J£ l(r, i) = ^ #f, i) 



and 

in terms of the bounded current and charge densities j and p, respectively. 



v W, *) + h & <K?> <) = 4 ^( f ' *) ( 4 ) 



It is convenient to define the fourier transform in time of a function f(r,t) by 

Ur) = ^Jdte^f(r,t), (5) 
so that wave-equations are transformed to the Helmholtz equations, with k = uj/c, 

and 

-(V 2 + k 2 )Mr) = ^ PuJ (r) ■ (6) 

We solve for the potentials in terms of the sources, and then reconstruct the electric and magnetic 
fields. 

The approach described here has some overlap with that of Rose [1]. 
An arbitrary vector field, A{r) can be expanded in spherical waves 

oo J 

A(f) = Y J E AW,?). (7) 
j=o M=-J 

The expansion coefficients can be given in terms of vector spherical harmonics. 

A(J, M, r) = i ]T [MJ, M; r)Yj&] (8) 
i=j-\ 

where fe(J, M; r) is a function of the radial coordinate r, and Yjf x are the vector spherical harmonics 
defined by 

?Wi(M) = E Y,{imlq\tUM)Yr{eA)e\. (9) 

m=—i q=— 1 

In (9) is a spherical unit vector, Yf l (9,<f>) is the usual spherical harmonic satisfying Yf 1 * = 
(— l) m Y™, and (£ml\£UM) is a Clebsch-Gordan coefficient. If e x , e^, e 2 , are the three rectangular 
unit vectors, then 



e\ — ^ (e x + ie y ) 



67 1 = 



e = e 
Q* _ 



(-l)«e_ s . (10) 
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The vector spherical harmonics obey the orthogonality 



dtt 



Y 



V M' 
' 1 Ji'l 



Sjj'Sw6mm' ■ 



Using this one finds 



± f t (J, M;r) = J dn [Y&(0, </>)] * • A(r) . 



(11) 
(12) 



In Lorentz gauge, the vector potential can be written in terms of the Green's function 

A,(r) = \ J d\'Ur')G k {\r- r'\) , k=^ (13) 

where G k {\r — r'\) = e '^J,\ ■ Hence, for (13) one has 



lf e (J,M;r) = ^ J Jdn [Y J i 1 (e,cl>)]*G k (\r-r 



(14) 



Consider the expansion of (14) for bounded sources, for which r' < R and r > R, 
j dSWji?G k {\r-r'\ 



,. £ 1 

dn Yl E (£mlq\eiJM)*Yr*(0,4>)eTG k (\r-r'\) 

J m=-t q=-l 



oo t 1 



53 (£mlq\eiJM)*eT f dnY E m *(.0,<t>)J2 E ikh$\kr)x 

i=-i q=-l J £'=0 m'=-V 



m=—£ q 

u(kr')Yp'(e,<t>)Yp'*(e',<i>') 
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53 53 (£mlq\eiJM)*el*Y e m {e' ^'Yiikh^ {kr)j e (kr')} 

m=—£ q=—l 

Yffi{#J)ikhP{kr)U{kr'). 



(15) 



In summary 



' i/*(J,M;r) = ikhi 1 \kr)lJd^r'j e (kr')j UJ (f')-Y 3 f 1 *(6', 
J+i 

AUJ,M,r) = I ^ f e (J,M;r)YjM(e,4>) 
i=j-\ 

oo J 

^(f) = E E ^(«/,M,f). 

J=0 M=-J 



(16) 



Putting this together, we can write 

oo J J+l 



4,(0 = E 53 53 ikh^(kr)af ei Yj^(e,, 
j=o m=-j e=j-i 



(17) 



where the coefficient 



a jei 



= \ J d 3 r'j e (kr') [? u (r ') • YjfW, cj>') 



(18) 



or 



Similarly, for the scalar potential 

Mr)= J d 3 r' Pu) (f')G k (\f-r'\) 

no t 



where 



Mf) = E E ikh\ L) (kr)bTYr(e,tf>) 
e=o m=-e 

bf = J d 3 r'p(f')j e (kr')Y e m *(0'A') ■ 



(19) 
(20) 
(21) 



The four coefficients b™ and are not independent. They are related by the conservation of 
charge, which is satisfied by computing the Lorentz condition, which will also ensure the gauge 
invariance of our results for the fields E and H. So one must satisfy V • — ikcj)^ = in this 
gauge. Before doing this, it is convenient to tabulate some useful formulae. 



For any <J>(r) 



V 
V 
V 

V x 

V x 

V x 



Hr)Yj M j(9, q 
Hr)Yj M J+1 (9,<f>) 



$(r)Y M 



J,J+i 



§{r)Y M 



J, J -i 



V[<S>(r)Y e m (6,<P)} = 



J+i 

2J+1 



A + J+ 2 

dr r 



®{r)Yj M j-i{e^)] = y/ 



2J+1 



d_ _|_ J+2 

dr r 



J-l 



dr 



$(r)Y e m (6,0) 
$(r)y/"(0» 



2.J+1 1 J,J 



Y 



M 



d_ 

dr 



J-l 



= i 



A _ I 

dr r 



Y 



M 



2J+1 1 J,J+1 



A + 1+1 

dr r 



(T + ^)*«i5?-i(fi.' 



2£+l 



(22) 



Now we can use the following properties for any spherical Bessel function zi(p). 



ze-i(p) + ze+i(p) 



21+1 



ze(p) 



d z e {p) 



dp 



[&*_i(p)-(* + l)z*+i(p)] 



(23) 
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From these properties, we find 

V-[z J+1 (kr)Yj M J+1 

V x 



—k 



J+i 

2J+1 



zj(kr)Yj M (e,<t>) 



-k^^zj(kr)Yj M (6,, 



z J+1 (kr)Yj M J+1 ] = ikzj{kr)J^Yl 



M 
J 



V X 

V x 



'zj(kr)Yjfj] = ik[-z J+1 (kr)^^fjf J+1 



Using these equations, we find that the Lorentz condition implies, writing aju as aj^, 

°e — l \ a e,e+iy w+i a e,e-i\J 2Z+1 j • 



The magnetic field is i7 w = V x A^. Direct computation shows that 

00 J 



it* = E E ^{^(^ 

J=0 M=-J 



M 
J, J 



2J+1 U J,J+1 



J+l _M 
2J+1 "J,J-1 



/ 27TT /l J+i(^ r )^j!S+i + hjl^kr)^ Yjfj + i\ } • 
The electric field is constructed from E u = —V<j) w + ikA^. 
Again by direct computation, and using (25), we have 



k = w 2 E E 



J=0 M=-J 



2J+1 U J,J-1 



hjl^kr)^ 27^]- Yjf J+1 - hjl^kr)^ Yjfj-i 



One can verify that all the Maxwell equations are satisfied. 

We can define the magnetic multipole \jJf = a^j, and the electric multipole 



Pf 



2.J+1 U J,J+1 



J+l n M 
2J+1 a J,J-l 



With these multipole coefficients, the fields take on a more compact, and symmetric appearance. 



oo J 
J=0M=-J 



M 
J 



E„ 



oo J 
J=0M=-J 

Note that if pf < — > nf, then 



(29) 



M 

,J+1 



i£ (electric) 
.E (magnetic) 



< — > 
< — > 



—H (magnetic) 
+H (electric) 



where the coefficients of p^f and pJf are denoted 
(electric) and (magnetic) respectively. 



One can put some of these terms in a more familiar form, if one notes that 

I J n M I J+l n M 

- y 2j+t a J,j+i ~ v 2j+T q j,j-i 

= \ J d 3 r'Ur') ■ { VOT n + i(kr')Yfff +l - ^ jj-^kr' 

= %f d 3 r'Un • V x [jj(kr')Yj M f} 

= a?/ ^bVMi^y] • [V x jL(r)] . 



(30) 



One can prove the identity 



£ ~ r 1 i ~ V W + V 2Z+T • 



(31) 



Using this, and the previously noted identities 
V x 



Hence 



jj(kr)fYj I 



J d 3 r'Ur') 



(32) 



V x 



jj(kr%Yj M *(9',cb') 



d 3 r'[V x (r'Ur '))] ■ [jj(£r')£ 1^(0', 



(33) 



^ | dV[V x • f' [jj(fcr')Y> M *(e',^) 



Note that only the transverse component of the current enters as the source, although in the 
longwavelength limit the longitudinal component of the current is the source for the electric 
multipole moments. This point is discussed in more detail by French and Shimamoto [3] and 
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by Snowdon [4]. Let us seek some alternate forms for our multipole expansion. Note that 
V x [zj(kr)rYj M ] = ^ J{J + 1) Zj {kr)Yfj, so that 



zj(kr)Yfi = ^fg=^ x [zj(kr)fYj 

zj(kr)Yj M 



M 



-L 



where L = —irx V. Also 
zj+i 



fVx 



zj(kr)Yj M j^ 

zj{kr)Yj M 



Using this 



that is 



Pf = % I d*r 



3 / 



-1 VxL 



zj(kr)Y^ 



p M 



\fjUW) 



jj(kr')Yj M * (9', #)} • [Vx jL(r') 



Vx JL(F') ■ L ^(Ax')^*^') 



and 



oo J 



oo J 

4 = w 2 E E 



(34) 



(35) 



(36) 



(37) 



J=0 M=-J 



pfijVxL) 
y/J(J+l) L 



4 1} (fcr)yf (^^)] + -^f= [^(^/(M)]} . 



(38) 



Our results are now expressed in the same form as French and Shimamoto [3], as expected from 
the gauge invariance of the electric and magnetic fields. 
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